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ABSTRACT 

7 -ray spectra of pulsars have been mostly studied in a phenomenological way, by 
fitting them to a cut-off power-law function. Here, we analyse a model where pulsed 
emission comes from synchro-curvature processes in a gap. We calculate the variation 
of kinetic energy of magnetospheric particles along the gap and the associated radiated 
spectra, considering an effective particle distribution. We fit the phase-averaged and 
phase-resolved Fermi-LAT (Large Area Telescope) spectra of the three brightest 7 - 
ray pulsars: Geminga, Crab, and Vela, and constrain the three free parameters we 
leave free in the model. Our best-fitting models well reproduce the observed data, 
apart from residuals above a few GeV in some cases, range for which the inverse 
Compton scattering likely becomes the dominant mechanism. In any case, the flat slope 
at low-energy (< GeV) seen by Fermi-LAT both in the phase-averaged and phase- 
resolved spectra of most pulsars, including the ones we studied, requires that most of 
the detected radiation below ^GeV is produced during the beginning of the particle 
trajectories, when radiation mostly come from the loss of perpendicular momentum. 


1 INTRODUCTION 

The wealth of F ermi-LAT (Large Area Telescope) data 
dAbdo et al.ll 201 ^ 'l has boosted our knowledge about 7 -ray 
emission from pulsars, allowing a better understanding of 
the fundamental high-energy processes which are responsi¬ 
ble for the conversion of the rotational energy into radiation. 
Two main channels are the likely origin of the detected radi¬ 
ation: photons emitted by particles moving in curved mag¬ 
netic fields and accelerating electric fields, i.e., the synchro¬ 
curvature (SC) radiation, and the inverse Compton (IC) 
scattering of backg round photons against energet ic magne¬ 
tospheric particles (IBogovalov fc AharoniarJl 200 Cfl . 

This paper is t he continuation of a series of works 
dVigano et al.l l2015al lalal. in which we focus on the high- 
energy SC radiation, and to which we refer the reader for 
further details. Our approach applies to any general gap 
located in the outer magnetosphere (even outside the light 
cylinder), since our parameters do not depend on an a-priori, 
detailed choice of the location of the gap. Instead, we show 
here how data can be used to constrain the values of the 
most relevant gap parameters of the model. 

The E^dN/dE spectra of pulsed 7 -ray emission from 
most pulsars peak around a few GeV, above which the hux 
quickly decreases with energy. Spectra are usually described 
by a cut-off power-law: 

^^PoA^exp 

where the four parameters to be fitted are the normalization 
Po, the low-energy slope index jj,, the peak energy Ep, and 
the exponential index s. The values of the best-htting models 



give a /i-range in the interval (—1,0.5), wi th the distribution 
peaking close to /r ~ —0.5 (see Fig. 7 of lAbdo et al.l[20131 . 
where their F is related to ^ as F = 1 — fj.). This relative 
flatness at low energies (0.1-1 GeV) contrasts with the value 
of /i predicted by the SC spectrum of a single-particle, which 
has a value /r = 0.25 due to mathematical properties of 
the functions involved. The latter is consistent only with a 
small minority of observed pulsars. This fact means that a 
simple SC radiation model which considers a mono-energetic 
distribution of particles is unable on a hrst-principle-basis to 
explain most of the observed spectra. 


The relative hatness (i.e., /r ~ — 1) of many pulsar spec¬ 
tra at low Permi-LAT energies is a basic issue that to our 
knowledge has never been properly addressed. In partic¬ 
ular, we believe it could signal that non-saturated parti¬ 
cles are important contributors to the total emitted radia¬ 
tion. Put otherwise, that the contribution of different parts 
of the particle trajectories is non-uniform. In IVigano et al.l 
(l2015all3j we showed how a large weight given to the ini¬ 
tial parts of the trajectories, where the radiation is domi¬ 
nated by synchrotron-like emission, can explain flatter slopes 
{fj. < 0.25). We study this possibility in more detail here, and 
advance that it is in agreement with data. 


Pulsar spectra show another important feature: for 
some sources, the phase-averaged spectrum shows a sub¬ 
exponential cut-off, i.e., s < 1. The high-energy tail (E > 
few GeV) decreases in energy slower than the expected SG 
radiation emitted by mono-energetic particles, which pro¬ 
duces a purely exponential cut-off, s = 1. This issue alone 
does not rule out SG radiation as the dominant mechanism: 
the phase-averaged spectrum is the superposition of radia- 
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tion coming from different phases, each of which could be 
properly described by a SC spectrum. The superposition of 
such spectra, with different energy peaks, slopes and fluxes 
can well lead to a sub-exponenti al cutoff. Such sce nario is 
found for the Geminga pulsar fsee lAbdo et al.|[2010al and the 
discussion below). The study of the phase-resolved spectra 
is the only way to assess this issue, and possibly rule out the 
SC radiation as the dominant mechanism above a few GeV. 

In this work, we calculate the variation of energy and 
the radiated spectra of magnetospheric particles along a 
gap, considering an effective particle distribution. We fit 
the phase-averaged and phase-resolved Fermi-LAT spectra 
of the three brightest 7 -ray pulsars: Geminga, Grab, and 
Vela. With such models we address the above mentioned 
issues. In Sj2] we recall the formulae describing the SG ra¬ 
diation, the computation of the particle trajectories and 
the associated radiative losses spectra. In |J3] we fit the 
phase-averaged and phase-resolved spectra of the radio-quiet 
PSR J0633-H746 (he reafter, Geminga as it is also known; 


Hal^ern^J^JHolt 1992 ), PSR J0534-I-2200 (the Grab pulsar; 


Staelin fc Reifenstein 196^), and PSR J0835-4510 (the Vela 
pulsar: Iharge et al.lll968fK By comparing models with data, 
we constrain the key physical parameters of our gap models, 
in particular, the component of the electric held parallel to 
the magnetic held (hereafter, i?||) and the ehective particle 
distribution. In S|4l we discuss the obtained results and its 
implications. 


2 SC RADIATION 

2.1 Single-particle radiation 

The SG radiation emitted by a single, charged particle spi¬ 
ralling around a magnetic held line is desc ribed by the fol¬ 
lowin g formula for the energy spectrum fsee lCheng fc Zhans 


a md the subsequent reformulation by IVigano et al. 
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In the formulae above, m and P are the rest mass and the 
Lorentz factor of the particle, a, rgyr and are the pitch 
angle, the Larmor radius, and the radius of curvature of 
its trajectory, respectively, e the elementary charge, B the 
local strength of the magnetic held, h the reduced Planck 
constant, c is the speed of light, Kn are the modihed Bessel 
functions of the second kind of index n, E is the photon 
energy. Pc is the characteristic energy of the emitted radia¬ 
tion. The factors Tefr, Q 2 , f, and Qr are introduced to pro¬ 
vide a compact formula. In the limits of high or vanishing 
perpendicular momentum, Eq. m reduces to purely syn¬ 
chrotron (^ 1 ) or curvature radiation (^ -C 1 ) formulae, 

respectively. In any case, the peak of the spectrum is located 
close to Pc and, for energies P ^ Pc, the do minant term in 
Eq. El , F{y), provides dPsc/dE ~ See IVigano et al.l 

(l2015ah for a more detailed discussion. 


2.2 SC radiation from the pulsar magnetosphere 


The detectable radiation coming from the magnetosphere 
can be estimated by integrating the single-particle spectrum, 
Eq. ([2]), along the travelled distance between the assumed 
boundaries Xin and Xout, convoluted with an effective par¬ 
ticle distribution, d N/dx, where x is th e distance along the 
magnetic field line dViganb et al.ll^lSd') : 


dPgap ^ r°“‘ dPsc dN^^ 

dE-y dpy dx 


( 12 ) 


The effective particle distribution represents the number of 
particles, per unit of distance, emitting radiation towards 
us. In general, it differs from the total particle distribution 
which is very hard to be inferred on purely theoretical basis. 

There exist two possible approaches to infer such a 
distribution. The first one is a sem i-analytical one, fol¬ 
lowed by e.g., I Zhang fc Chengl lll997ll . in which most gap 
parameters are univocally determ ined by the values of P 
and P*. In IVigano et al.l (l2015bll , we questioned such ap¬ 
proach, analysing the impact of its many underlying as¬ 
sumption and approximations, which necessarily relies on, 
for instance, taking single-value parameters in many quanti¬ 
ties, neglecting their uncertainties and dispersions (caused, 
e.g., by viewing angle, beaming effects and the magneto- 
spheri c geometry), or fixing the pitch angle to a non-zero 
value l|Zhang fc Cheng|l997| ), whil e it should go to zero un¬ 
der SG losses JVigano et al.lf^l5al ~l. 

On the other ha nd, the numerical s i mulat ions by 
iTakata et al.l (l2006l j and lHirotanil (l2006l . [20071 . l20ni^ 'l include 
the pair cascade, the Boltzmann equation for particles and 
the interaction with the radiation. These simulations are 
physically more self-consistent, but they also rely on some 
intrinsic, hardly avoidable assumptions, like the precise loca¬ 
tion of the gap, the intensity of incoming/out going currents 
(if any), and a simplified magnetic field geometry. Moreover, 
the computational cost makes this approach unviable if one 
wants to fit real data once a nd again. More recent state- 
of-the-art MHD simu lations l|KalaDotharakos et al.l l2014l : 
iBrambilla et al.l l2015l j roughly predict the curvature radi¬ 
ation emitted outside the light cylinder, under certain as¬ 
sumptions (negligible pitch angle, spin period of ~ ms, con¬ 
stant magnetospheric conductivity outside the light cylin¬ 
der, and purely force-free magnetosphere inside it). 
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Figure 1. Evolution of the Lorentz factor F (top), the SC param¬ 
eter ^ (middle), and the single-particle SC power Psc (bottom) 
calculated for a model with the timing properties of Geminga 
pulsar (see Table (2]), and varying one by one the relevant param¬ 
eters, as indicated in the legend, over a baseline model consisting 
of: _E|| = 10^ ®*^ V rf = 0.5, b = 2.5, Fin = 10^, dm = 7r/4, 

^\n/R\c = 0.5. Some of the lines overlap with the base model 
(black solid line). 

These difficulties in fitting spectra from first principles 
led us to take a more effective approach, based on a gap 
model containing a few effective parameters, which allow a 
systematic study of 7 -ray spectra from which to extract few, 
but solid, physical conclusions. We explore here the effects 
of having effective particles distributions with larger weights 
in the inner par t of the trajec t ory (sm all xq/R\c: see below), 
as discussed in IVigano et al.l (|2015al l3). This choice reflects 
both an observational and a theoretical consideration. As 


Table 1. Parameters entering in the SC model, explored range, 
and treatment in the fitting procedure. For a given pulsar, P and 
P* are fixed by timing observed properties. The second group 
of three parameters are the ones constrained by data (we indi¬ 
cate within parenthesis which is the most constraining spectral 
property), while the remaining six parameters have much smaller 
influence on spectra and are fixed to the listed values. 


Parameter 

Range 

Treatment 

p 

- 

timing data 

B„ 

- 

fixed by timing data 

log(A||[V m-l]) 

6.5-9.5 

fit (energy peak) 

No 

1026_io34 

fit (luminosity) 

xo/Ric 

0.001-1“ 

fit (low-energy slope) 

V 

0.2-1.0 

fixed to 0.5 

b 

2-3 

fixed to 2.5 

^in/^lc 

0.2-1.0 

fixed to 0.5 

3^out / R\c 

1.0-2.0 

fixed to 1.5 

Tin 

lO^-lO"' 

fixed to 10® 

Oin 

0 

4 

to 

fixed to 7r/4 


“ See Eq. (I13II for the definition of the effective parameter distri¬ 
bution dN/dx. We also explore the uniform particle distribution, 
dN/dx = No/(xout — ^in)- 


said in the introduction (see also lVigano et ai]|2015(Jl a uni¬ 
form distribution is not compatible with having the variety 
of slopes at low energies, E < GeV, seen by Fermi-LAT. 
Physically, we can consider a number of effects that can 
qualitatively justify a large weight in the low ^ parts of the 
trajectories: the beaming is larger at smaller P, so that it 
is easier to detect radiation for less energetic particles. The 
geometry of the emitting region can also affect the radia¬ 
tion. The cascade of pairs, produced especially in the inner 
gap, where the interactions with X-ray photons are more 
likely, can also give a net effect of enhancing visible emis¬ 
sion from particles with ^ 1 , when they lose perpendicular 

momentum. Other effects, related to magnetic geometry and 
viewing angle, can further complicate the picture, and make 
any theoretical prediction quite difficult. 

In order to introduce such effective particle distribution 
in a way that allow fitting, we shall introduce a functional 
form depending on the length-scale xo and normalization No 
(total effective number of particles): 

dN e-(^-^in)No 

—— = No — 7- -- r . (13) 

dx * 0(1 — 

The choice of this particular function is motivated by sim¬ 
plicity. Other functions, with other parameters allowing to 
produce a non-uniform particle distribution along the gap 
could have been considered. We stress that this choice is 
an effective way to parametrize our ignorance about the 
combined effects described above on the observed emission, 
and will allow to constrain significant qualitative features 
emerging from the data fits. We also stress that Eq. m 
shall represent the particles which effectively produce the 
observed radiation, not the whole particle population in the 
gap, which we may never infer. 
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Table 2. Pulsars considered in this work: timing properties, references for the spectral analysis, phases considered (average, peaks and 
minimum or inter-pulse state), and the values of the relevant parameters of our SC best-fiting models, with the associated value of x^. 
The loss of the rotational energy is calculated as ^rot = ^ttIP/P^ = 3.9 x 10‘^®P/P[s]^ erg/s, while the characteristic age is defined as 
r = P/2P. 


Pulsar 

P 

P 

P* 

T 

Prot 

Refs. 

Phases 

log(P||) 

3^0 

No 

X® 


[s] 

[10-14] 

[lOi^G] 

[yr] 

[erg/s] 



[V m.“l] 

[Pic] 

[10®°] 


J0633-I-1746 

0.237 

1.1 

3.3 

1.2 X 10® 

3.2 X 10®4 

[1] 

average 

7.65 

0.013 

18.8 

18.1 

(Geminga) 







PI: 0.637-0.643 

7.80 

0.013 

0.26 

0.68 








P2: 0.131-0.141 

7.65 

0.010 

0.82 

0.91 








Min.: 0.942-1.000 

7.20 

0.026 

0.76 

1.17 

J0534-I-2200 

0.033 

42.3 

7.6 

1.1 X lO"! 

4.6 X 10®® 

[2] 

average 

8.75 

0.003 

41.3 

4.55 

(Crab) 







PI: 0.987-0.993 

8.60 

0.004 

3.16 

1.48 








P2: 0.366-0.386 

8.80 

0.003 

1.8 

1.48 








Bridge: 0.098-0.286 

8.45 

0.010 

2.9 

0.54 

J0835-4510 

0.089 

12.5 

6.8 

3.4 X 10® 

6.9 X 10®® 

[3] 

average 

8.05 

0.010 

14.6 

136 

(Vela) 







PI: 0.133-0.135 

8.05 

0.008 

0.22 

2.12 








P2: 0.562-0.563 

8.00 

0.014 

0.13 

2.76 








P3: 0.315-0.324 

8.15 

0.013 

0.076 

1.39 


References: [ll lAbdo et al.l l|2010ah : [2l lAbdo et al.l ll2010bl') : [3l lAbdo et al.l l|2010d) . 


2.3 Particles trajectories 

The quantities defined in Eqs. dl-lIB, appearing in Eq. (H, 
depend on the local values of F, a, r^, and B. The kinematic 
parameters, F and sin a, evolve along the gap due to elec¬ 
tric acceleration and radiation losses, thus we simulate the 
motion of charged particles by evolving the parallel and per¬ 
pendicular momenta_ofpajticles according to the equations 
of motion JVigano et al.llioiSal h 

d(psina) PacSina 

dt V 

dipcosa) ^ Psccosa 

-;- = ePii- 

dt " V 

where p = Tmv is the momentum of the particle, v its 
spatial velocity, i5|| is the parallel componen t of the elec¬ 
tric fie ld, considered constant along the gap dVigano et alJ 
l2015bl l. and Psc is the single-particle SC power, obtained by 
integrating Eq. © in energy: 


(14) 

(15) 


2e^F^c 

3r2 


ffr • 


(16) 


On the other hand, we prescribe simple functional forms of 
rc(x) and B(x), de pending on two effect ive parameters, rj 
and h, as justified in I Vigand et al.l ll2015bl l3l: 

(u) 

B{x) = P* , (18) 

where P* is the neutron star (NS) radius, and we employ the 
standard definitions of light cylinder distance and inferred 
dipolar magnetic field at the polar surface: 

Pic = 4.77 x 10®P[s] cm , (19) 

27r 

P* = 6.4 X 10^® y^P[s]P G . (20) 


Expected valu es for a pulsar magne tosphere are rj ~ 0.2 — 1 
and b ~ 2 — 3 (IVigand et al.l[^15lJ) . 

We summarize he re the main results obtained in 
our previous works dVigano et al.l I2015al l3. see also 
iHirotani fc Shibatall999h in Fig.[T] by showing the evolution 
of the Lorentz factor F (top), the SC parameter ^ (middle), 
and the single-particle SC power Psc (bottom), obtained for 
models having the timing parameters of the Geminga pul¬ 
sar. Soon after pair creation, ^ S> 1, i.e. the losses mostly 
regard the perpendicular momentum, being approximated 
by synchrotron formula. Then, in a length-scale x -C Pic, 
the electrical acceleration makes F increase and sin a de¬ 
crease by orders of magnitude, until the radiative losses can 
be approximated by purely curvature radiation, ^ <C 1. In 
this regime, F and Psc tends to saturate to a value deter¬ 
mined by the balance between the electric force and the 
radiation-reaction force. 

The baseline model (black line) has the following pa¬ 
rameters: P = 0.237 s, P* = 3.3 X 10^^ G (as in the Geminga 
pulsar), P|| = 10^ ® V m~^, 77 = 0.5, b = 2.5, im = O.SPic, 
Fin = 10®, and Qin = 7 r/ 4 . The other lines show models ob¬ 
tained exploring reasonable variati ons of the relevant par am- 
eters, according to what studied in I Vigand et al.l (l2015bll . 6 , 
Fin, and min have visible effects only in the very early part of 
the trajectories, x/Ric < 10“®. This happens because, when 
5 < 1 , then the radiation is approximated by curvature ra¬ 
diation, which is independent on the value of P and a. The 
parameters Xin and rj have an overall negligible influence as 
well, limited to minor differences in the outer part of the 
trajectory. Hence, the most relevant parameter affecting the 
trajectories is, by far, Py: the larger it is, the larger F and 
the smaller 

In conclusion, we can see that for x/Ric < 10“® — 10“ ® 
(value depending on the model) the deviations from a purely 
curvature radiation (^ <C 1) can be important. This is why 
we need to consider the full expression of SC radiation and 
evaluate it along the whole trajectory of the particles mov- 

















































ing through the gap, especially if most of the detectable 
radiation originates at the beginning of their trajectories, 
when particles are relatively slow and radiate mostly due to 
perpendicular momentum losses. 


3 RESULTS 

Instead of the (sub-)exponential cut-off power-law fit, 
Eq. 0, we aim at fitting the spectra with the expected SC 
radiation, constraining the physical parameters of our gap 
model. In order to quantify the relative goodness of the fits, 
we compare the observed Eermi-LAT data with the binned 
theoretical spectra. In Appendix|0we give the details about 
the explored grid of parameters, the re-binning of spectrum, 
and the definition of the goodness-of-£t indicator x^. 

According to the model setup described above, for a 
given pulsar we can fix the values of the spin period P, 
which dehnes Ric, Eq. m, and the surface magnetic held 
R*, inferred from Eq. (I 20 j. Then, we have six parameters 
that regulate the trajectory: r), b, Xin, Tin, and ain, and 

three parameters describing the effective particle distribu¬ 
tion, Eq. 0 : No, xo/Ric and Xout/Ric 

Since our spectral models are not analytical, and re¬ 
quire the computation of the particle trajectory, we cannot 
perform a systematic, large and blind coverage of the mul¬ 
tidimensional space of parameters. Instead, as a hrst step, 
we start by identifying which parameters have the largest 
impact on the spect rum and focus on th ose, relying also on 
our previous works dVigano et al.|[^15al l3i. 

In Fig. m we saw that the kinematic evolution of parti¬ 
cles is very sensitive on the value of i?||, while the other Hve 
parameters, {ri, b, Xin, Fin, ain), have non-negligible impact 
only in the very limited parts of the trajectories. In partic¬ 
ular, b, Fin, and ain will have an impact only if the particle 
distribution has a width xo/Ric < 10 “'*. 

Thus, since the typical photon energy, Eq. ®, strongly 
depends on the value of F, E|| is mostly constrained by the 
energy peaks of observed spectra. The value of xo/Ric is 
mostly constrained by the low-energy slope. The value of 
Sout has an effect only if xo Sout, for which the outer 
part of the gap give important radiation. However, we an¬ 
ticipate here that small values xo/R\c 1 are favored from 
the spectral hts, thus Xout is unconstrained. Lastly, the nor¬ 
malization No is picked up by minimizing x^, once the other 
parameters are Hxed (see Appendix [0. 

Briehy: the timing-derived parameters, P and R*, are 
Hxed for a given pulsar, and the three most relevant param¬ 
eters are: No, xo/R\c and Ry. In Table [T] we summarize the 
set of parameters and the treatment we use when we Ht data. 
If not specihed, the other parameters are Hxed to a baseline 
model with the values indicated in the third column. 

Below, we apply our models to Ht data to Geminga, 
Crab, and Vela pulsars, considering both the phase-averaged 
and phase-resolved spectra of Rermi-LAT. In Table 0 we 
show their timing properties, the derived rotational energy 
loss, characteristic age and surface magnetic Held, the phases 
for which we consider the spectra below, appropriate refer¬ 
ences, and the results of the Ht that we illustrate below. 
Since we focus on the SC radiation, we limit our study to 
the Eermi-LAT data. 
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3.1 Geminga 

The Geminga pulsar has been detected in 7 -rays and X- 
rays. Upper limits on the E > 100 GeV emissi on have re¬ 
centl y ruled out any hardening above > 50 GeV (lAliu et al.l 
l2015li . The 7 -ray light curve of the Geminga pulsar is double- 
peaked, with PI being the main peak for E > 0.3 GeV. The 
second peak, P2, is slightly higher than PI in the range 0.1- 
0.3 GeV. For the phase-resolved spectra, we consider the 
bins associated with the two peaks and the mini mum phase, 
identiH ed according to the Huxes reported by lAbdo et al.l 
(l2010al ). and indicated in Table [21 

We start the discussion with the phase-averaged spec¬ 
trum. In the top left-hand panel of Fig. [2] we show the 
best-Htting model and the observed spectrum, obtained by 
exploring the log(E|| )-log(a:o/Ric) plane, Htting No in each 
model, and Hxing the other parameters to the baseline model 
of Table[T](see above and Appendix|0for more details about 
the Htting procedure and grid of explored parameters). The 
best-Htting values are Ey = 10^'®® V m“^, xo/R\c ~ 0.013, 
and No = 1.9 x 10®^ (see also Table[21). In the top right-hand 
panel, we show the related contour plots of X^/Xmim where 
Xmin is tile minimum best-Htting value (white cross). A low 
value of xo/Ric is needed to Ht the spectrum below the en¬ 
ergy peak, which require that more weight is given to the 
parts with large i.e., dominated by synchrotron-like emis¬ 
sion. This is a clear example of why a uniform distribution, 
where the radiation would be dominated by the ^ -C 1 part 
of the trajectory, cannot explain the low-energy slope. 

The Xinin associated with the best-Htting of the average 
spectrum is large for three reasons. First, the errors in the 
low energy part, < 1%, are very small. Secondly, the high- 
energy (above 10 GeV) data points are not Htted well. This 
incompatibility is inherent to the strictly exponential decay 
of the SG radiation formula, and agrees with the phenomeno¬ 
logical sub-exponential cut-off model claimed in lAbdo et al.l 
(l2010al ). We come back to this below. Thirdly, further Hne- 
tuning of the explored parameter (e.g., a Hner grid of Ey 
and xo/Ric) would further improve the Ht. 

Looking at the contour plot, we note a strong anticorre¬ 
lation between Ey and xo/Ric, with a strip of models having 
yf < 2xf,in that constrain log(xo/Ric) and log(Ey) within 
~ 20 and ~ 5 per cent, respectively. This anticorrelation, 
which is seen in all cases we have studied, is interpreted as 
follows. The low-energy slope with < 0.25 can be repro¬ 
duced only if the initial parts of the trajectory are domi¬ 
nating the observed radiation. A large contribution of such 
parts, 5 ^ 1 , can be achieved by different models: the larger 
the electric Held, the shorter is the acceleration length-scale, 
hence a smaller value of xo/R\c is needed to consider the 
parts with relatively low F. In any case, the favored values 
are xo/Ric 1: we are seeing radiation only from the inner, 
initial part of the particle trajectories. 

As a second step, we check the inHuence of other param¬ 
eters. We evaluate the best-Htting spectra for different val¬ 
ues of xo/Ric, with Ey = 10^'®® V m“^ (best-Htting value), 
and varying, one by one, the following three parameters: 
b £ [2,3] (step 0.1, left-hand panel), 7 £ [0.2,1.0] (step 0.1, 
middle), Xin/Ric £ [0.2,1.0] (step 0.1, right-hand panel). Re¬ 
sults, drawn in the bottom panels of Fig. (21 show that these 
three parameters are hardly constrainable, apart from ap¬ 
parently favoring Xin > 0.5Ric and b > 2.5. Both weak con- 
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Figure 2. Study of the phase-averaged Fervnj-LAT spectrum of the Geminga pulsar. Top left: Fermi-LAT data and corresponding 
best-fitting model (see also Table[2]l. Top right: contour plots for the values of X^/Xmin ™ log(F||)-log(a;o/Fic) plane, with the values 
of the remaining parameters given by the baseline gap model, see Table [T] The white cross represents the best-fitting model of the 
left-hand panel. Bottom panels: similar contour plots for the planes 6-log(a;o/iJi;;)i 7;-log(a;o/Fic), and 3;in-log(a;o/Fic), respectively, with 
log(F|| [V m~^]) = 7.65. See text and Appendixfor details about fitting procedure. 


strains correspond to models with more efficient radiative 
losses: low values of b correspond to larger values of B in 
the outer magnetosphere, low values of Xin imply a smaller 
radius of curvature. In both cases, particles will reach lower 
values of F. However, the large degeneracy and the much 
stronger dependence on i7|| and Xo/Ric prevents us from 
taking firm conclusions about these lower limits. The lack 
of such constraints justifies our initial choice of iJ||, xo/Ric, 
and No as the only three parameters to be left free in the 
fit. 

Since phase-resolved spectra show important variation 
between phases, we explore three representatives of them 
here. In Fig. [3] we show the same best-fitting model and 
contour plots shown for the phase-averaged spectrum, but 
referred to the spectrum of PI. In these cases, all data can be 
fitted well by our models, giving a much smaller (mainly 
due to the fewer photons, i.e., larger errors), and fitting well 
also the spectrum above GeV. The correlations between pa¬ 
rameters are similar to the ones of the phase-averaged spec¬ 
trum, with the only differences in the weak constraints that 
can be put on the less relevant parameters, Xi^, rj and b. In 
Fig. |4]we show the best-fitting models for P2 and minimum 
phases, which are compatible with observations. 


In Table [2] we list the best-fitting parameter values for 
all cases. Note that in the phases of the peaks the fitted val¬ 
ues of i5|| and xo / R\c are similar to those obtained with the 
phase-averaged spectrum, since the peaks are the main con¬ 
tributors to the average. On the other hand, the minimum 
shows a much softer spectrum, which can be reproduced by 
smaller i5|| and larger xo/Ric- 


In general, we conclude that the high-energy phase- 
resolved spectra of Geminga are compatible with SC models, 
if the radiation mostly comes from the loss of perpendicular 
momentum at the beginning of the trajectories. This agrees 
with the observational conclusions c oming from the deta iled 
phase-resolved analysis reported bv lAbdo et al.l ll2010ah : at 
each individual phase, the spectrum can be reasonably de¬ 
scribed by a simple cut-off power-law, with energy peaks, 
slope indices, and fluxes varying by a factor of a few between 
different phases. We confirm that the sub-exponential cut¬ 
off seen in the phase-averaged spectrum is given by the effect 
of summing up the spectra coming from different phases. 
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Figure 3. Same as Fig. [2] but for the phase-resolved spectrum in PI of the Geminga pulsar. 


Geminga P2 


Geminga Minimum 




Figure 4. Best-fitting models of the Geminga pulsar for P2 (left) and minimum (right). 


3.2 Crab 


The Crab pulsar has been detected in all energy ranges, 
from radio up to v ery high energies. Besides Fermi-LAT 
dAbdo et al.ll2010bl ). pulsed detections have been reported 


up to 2 TeV by MAGIcQ llAliu et al.l l2008l: AdeksidetaJj 
|2011| . [20 i3 i. above 100 GeV by VERITAS ( Aliu et al ■II2011I) . 
At all energy ranges, the light curve presents two well- 
defined peaks, PI and P2 of Table [2] with almost no lag 
between different phases. In our analysis, we consider the 


1 


http: //f ermi. gsfc.nasa.gov/science/mtgs/syniposia/2014/abstracts/185 
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Figure 5. Same as top panels of Fig. [2] but for the Crab pulsar. 

Crab P2 


Crab Bridge 





Energy [eV] 

Figure 6. Best-fitting models of the phase-resolved spectra of the Crab pulsar in PI (left), P2 (center), and bridge (right). 


phase-averaged spectrum, and the phase-resolved spectra 
for three phase bins, representing PI, P2, and the bridge 
between them. 

As before, we have noted that relatively tight con¬ 
straints can be put only on if||, xq/Ric, and Nq. This is 
a general feature, so hereafter we will neglect the other pa¬ 
rameters of the model. In Fig. [5] we plot the same contour 
plot and corresponding best-fitting spectrum to the phase- 
average data, as we did for the Geminga pulsar. Fig. [3 We 
find a satisfactory fit, with a constrained value of £1|| which 
has to be larger than that in Geminga, in order to reach the 
observed peak at few GeV. The value of xq /i?ic is very small, 
again indicating an important contribution by synchrotron¬ 
like emission. 

In Fig.[6]we show the best-fitting spectra for PI, P2 and 
the bridge. The best-fitting values (see Table o are similar 
to the phase-averaged spectrum, while the bridge spectrum 
requires a slightly larger xo /R\c and lower , because of the 
steeper low-energy slope. In general, the fit to phase-resolved 
spectra is satisfactory. 

Although we are able to reproduce the Ferroi-LAT spec¬ 
tra, we note that the recent detection of the pulsed emission 
in TeV (not considered here) rules out the SC model to 
be the only responsible mechanism, indicating a likely im¬ 


portant contribution of the IC mechanism at energies > 10 
GeV. 


3.3 Vela 

The Vela pulsar has been detected from radio to very high 
energies. Pulsed em ission has been detected above 50 GeV 
dLeung et al.ll20l3 . recently also by the H.E.S.S. telescopcQ), 
similarly to the Crab pulsar. It presents a large variability 
of light curves depending on the energy. In order to study 
the phase-resolved y-ray spectra, we consider data coming 
from the three phases listed in Table [21 corresponding to 
the two peaks, PI and P2, which are the maxima below and 
above 1 GeV, respectively, and an inter-pulse phase, P3, 
corresponding to a third peak visible at F > 3 GeV. 

In Fig. [7]we show the contour plot of X^/Xmim qualita¬ 
tively similar to the previous cases, and the corresponding 
best-fitting spectrum for the phase-averaged spectrum and 
its best fit for the Vela pulsar. As before, the tail at hard 
energy cannot be reproduced by SC models. The best-fitting 


^ http://www.mpg.de/8287998/velar-pulsar 
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Figure 7. Same as top panels of Fig. but for the Vela pulsar. 
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Figure 8. Best-fitting models of the phase-resolved spectra of the Vela pulsar in PI (left), P2 (center) and P3 (right). 


value of i5|| is intermediate between the Crab and Geminga 
cases, and again we need Xo/Ric <C 1 . 

In Fig. [ 8 ] we show the phase-resolved spectra for the 
three phases, and the best-fitting obtained within our mod- 
elsI 3 The spectra in PI and P3 are compatible with the SC 
emission model, while in P2 the residuals at high energies 
could point to contributions from IC scattering. 


4 CONCLUSIONS 

In this work, we have applied a SC model to the phase- 
averaged and phase-resolved Fermi-hAT spectra of the three 
brightest and best-studied 7 -ray pulsars: Geminga, Crab, 
and Vela. The theoretical spectra rely on a ID gap model 
based on effective parameters, whose exploration allow us 
to systematically £t data. For a given pulsar, P and 
are fixed by the observed timing properties. These define 
Ric, Eq. (O, and the surface magnetic field B*, Eq. (EOl)- 

^ The valu e s are n ot the same as reported in the Fig. 10 of 
lAbdo et al.l ll2010(J) . in which the flux values are renormalized 
by the inverse of the phase bin width (T. Johnson, private com¬ 
munication). 


Then, the properties of the emitting region can be effec¬ 
tively described by the following effective parameters: the 
accelerating electric field By, four geometrical parameters, 
77 , b, Xin/Ric, and Xout/Bic, two kinematic parameters. Fin 
and ctin, and two parameters describing the effective particle 
distribution: Nq and xo/Ric- 

A main conclusion is that By , xo /Bic, and No have much 
more influence on spectra than the remaining six parame¬ 
ters. Thus, observational data can quantitatively htted by 
our model varying just these three parameters (which are, 
in number, one parameter less than the sub-exponential cut¬ 
off power-law employed by Bermi-LAT to describe the data), 
keeping the others hxed. By is mainly constrained by the en¬ 
ergy of the peak, xo/R\c by the low-energy slope, and Nq by 
the flux. 

In Fig.[9]we show the evolution of the Lorentz factor for 
the best-fitting models of the three phase-averaged pulsars. 
All curves of r(a;/Bic) versus normalized distance along the 
line are similar, reaching similar values of F ~ few 10 ^ at 
a distance x ~ 10 “^ — a length-scale comparable 

with the inferred values of xq. 

In all cases, the best-fitting models are consistent with 
radiation mostly coming from the initial part of the particle 
trajectories (xo/Ric ^ 10 ”^), where the perpendicular mo- 
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Figure 9. Lorentz factors versus normalized distance along the 
gap, x/Ric, for the best-fitting models of the phase-averaged spec¬ 
tra of the three pulsars. 


mentum losses (i.e., synchrotron losses) are not negligible 
> 1, see Fig. [1]). This constraint is related to the rela¬ 
tively flat slope of the low-energy part, common to many 
other pulsars, and represents strong, qualitative indications 
about physical processes. In the framework of a SC interpre¬ 
tation of the pulsed emission, this can be related to several 
possible physical effects. The beaming angle of the emit¬ 
ted radiation is inversely proportional to F, thus low-energy 
photons, mostly coming from low-F particles, are easier to 
be detected. Moreover, the cascade process, not simulated in 
our effective approach, could actually provide a large num¬ 
ber of particles quickly losing perpendicular momentum in 
synchrotron-like radiation. 

For Geminga and Vela pulsars, our best fits to the 
phase-averaged spectra show residuals after the ~ GeV 
peak. However, most of the studied phase-resolved spectra, 
at each individual phase, are compatible with purely SC ra¬ 
diation. In particular, the resulting sub-exponential cutoff 
in the phase-averaged spectrum of Geminga is coming from 
the conjunction of individual exponentially decaying spec¬ 
tra showing different cutoff energies. As a consequence, it 
is natural not to expect any detection from Geminga above 
tens of Ge V, in agreem ent with VERITAS and MAGIG ob¬ 
servations (lAliu et al.ll2015|l . 

In some other cases, such as in at least one of the 
phases of the Vela pulsar, a sub-exponential cut-off results 
in important residuals at high energies even in the phase- 
resolved measurements. These results are consist ent with the 
report ed phenomenological fits by Fermi-LAT llAbdo et al.l 
|2010blHHl . Another qualitative difference among the three 
pulsars considered is that Vela and Grab pulsars show pulsed 
emission at several tens of GeV, and up to 2 TeV, respec¬ 
tively. Detections at such high energie s need to be explained 
by IC scattering. iLeung et al.l (l2014l f proposed a spectral 
model for Vela based on curvature radiation, consisting in 
the superposition of different outer gap structures, leading 
to different energy peaks. However, their model fails to re¬ 
produce the low-energy slope, and would be increasingly in 
tension the higher the energy at which pulsations are de¬ 
tected. 

Note that even such mismatch observed at large en¬ 


ergies in some cases does not invalidate the SG as the 
main mechanism of radiation below ~GeV. Instead, it likely 
points to extra-contributions coming from IG, dominating at 
large energies o nly in a few case s. Recent stacked searches 
for pulsations bv lMcGannl (|2014| ~) indicate no strong average 
emission at I? > 50 GeV, which could make the SG model 
valid for most pulsars, with Grab and Vela pulsars being ex¬ 
ceptions. Finally, note that the three pulsars span a range of 
two orders of magnitude in the estimated age (see Table [5)). 
The rotational energy differences could explain the different 
importance of the IG mechanism, which could be more ef¬ 
fective for young, more energetic and hotter NSs, due to the 
larger X-ray flux from the surface. 

Gomparing the obtained best-fitting values for the three 
pulsars, we note that the Grab pulsar needs larger values of 
ifII and No, compared with Vela and Geminga, for which the 
best-fitting values are similar. The similar values of No be¬ 
tween Vela and Geminga, in particular, are in contrast with 
the large difference (a factor ~ 200) between their rotational 
energies, and reflects the high efficiency of Geminga. 

A complementary study to our effective approach 
should ideally include 3D simulations of particle dynamics, 
considering the interaction with the radiation, leading also 
to the IG contribution estimate, with all the geometrical and 
beaming effects considered (and probably fitting on the gap 
location). Moreover, the inclusion of multiwavelength data 
and the study of the light curves should better constrain the 
models. 

In a forthcoming work, we will apply our models to 
systematically study the entire population of y-ray pulsars 
having good-quality data. Thus, we will constrain the values 
of if||, Xo/Ric, and No, and look for trends in the population. 
Since good-quality phase-resolved spectra are available only 
for a minority of pulsars, one can work with the best phase- 
averaged spectra, supported by the fact that, for the three 
cases here studied here and for which phase-resolved analy¬ 
sis exist, the obtained best-fitting values of i?|| and xo/Ric 
are similar to the phase-resolved spectra at the peaks (see 
Table El). 
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X = 


Abin — 1 


E 


/ T bin r bin \2 
\-^obs -^gapj 




(A2) 


where JLobs = is the luminosity error for each 

bin. Such error neglects the (likely large) uncertainties on 
distance and beaming factor, but this would only change the 
luminosity, thus providing a different inferred normalization 
No, but without changing the spectral shape, and therefore, 
the constrains on E|| and xo/Ric- 

In the plots of the paper, following the literature, we 
plot the binned functions E^dN/dE = in units 

erg s“^, both for the theoretical models and data. 


APPENDIX A: FITTING PROCEDURE 


The aim of the fitting procedure is to constrain the three 
relevant parameters of our model, E||, xo/Ric, Aq, by com¬ 
paring the theoretical and observational spectra, the lat¬ 
ter taken from the publicly available data of the second 
Fermi-LAT pulsar cataloglj Data consist of a number of 
Abin ~ 5 — 10 energy bins (each one with its extremes Ei 
and E 2 , and its weighted central energy Acent), each one with 
its associated photon flux, Ef^^ (in units photons cm“^s“^), 
and its associated statistical error, We consider the 

isotropic luminosity, L^b" = And^Ef^^, where d is the dis¬ 
tance of the pulsar to the Earth, and we neglect possible 
beaming effects which would reduce the inferred luminosity. 
We neglect bins having only upper limits to the flux. 

To explore the space of parameters, we span a grid of 
41 different values of A||, logarithmically equi-spaced, cov¬ 
ering two orders of magnitude, and a grid of 50 values of 
xo/Ric £ [0.001, 1] (with increasing steps), plus the case of 
uniform effective distribution of particle. For each pair of val¬ 
ues (A||, xo), we evaluate the expected spectrum over a grid 
of hundreds of points in the 100 MeV-100 GeV range. Then, 
we renormalize it by the best-fitting value of No, found by 
means of scanning a grid, having progressively finer steps, 
up to dNo/No < 10“^. Therefore, we integrate the SC pho¬ 
ton spectrum (related to the SC energy spectrum. Eg. 1121 by 
dAga.p/dA = (l/A)dPga,p/dA) in each bin (and normalizing 
by the bin width), to obtain the binned SC photon spectrum 
(number of photons per unit energy): 


L 


bin 

gap 


A2—Al 


1 dAgap 
A dE 


dA . 
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Therefore, we calculate the goodness-of-fit indicator in the 
standard way: 
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http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog/ 






